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Abstract 

Complex Hadamard matrices, consisting of unimodular entries with 
arbitrary phases, play an important role in the theory of quantum in- 
formation. We review bàsic properties of complex Hadamard matrices 
and present a catalogue of inequivalent cases known for dimension N = 
2, ... ,16. In particular, we explicitly write down some families of complex 
Hadamard matrices for N = 12, 14 and 16, which we could not find in the 
existing literature. 

1 Introduction 

In a 1867 paper on 'simultaneous sign-successions, tessellated pavements in two 
or more colors, and ornamental tile-work' Sylvester used self-reciprocial matri- 
ces, defined as a square array of elements of which each is proportional to its first 
minor pP . This wide class of matrices includes in particular these with orthogo- 
nal rows and columns. In 1893 Hadamard proved that such matrices attain the 
largest value of the determinant among all matrices with entries bounded by 
unity 0. After this paper the matrices with entries equal to ±1 and mutually 
orthogonal rows and columns were called Hadamard. 

Originally there was an interest in real Hadamard matrices, Hij € R, which 
found diverse mathematical applications, in particular in error correction and 
coding theory [3] and in the design of statistical experiments 0]. Hadamard 
proved that such matrices may exist only for N = 2 or for size N being a 
múltiple of 4 and conjectured that they exist for all such N. A huge collcction 
of real Hadamard matrices for small vàlues of N is known (see e.g. but 
the original Hadamard conjecture remains unproven. After a recent discovery 
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of N = 428 real Hadamard matrix [7], the case N = 668 is the smallest order, 
for which the existence problem remains operi. 

Real Hadamard matrices may be generalized in various ways. Butson in- 
troduced a set H(q, N) of Hadamard matrices of order N, the entries of which 
are q-th roots of unity [311] ■ Thus H(2, N) represents real Hadamard matrices, 
while -ff (4, TV) denotes Hadamard matrices 1 with entries ±1 or ±i. If p is prime, 
then H(p, N) can exist if N = mp with an integer m [H], and it is conjectured 
that they exist for all such cases [T4] . 

In the simplest case m — 1 Hadamard matrices H(N,N) exist for any di- 
mension. An explicit construction based on Fourier matrices 

[F' N ] k := -L e »ü-i)(*-D^ with j,ke{l,...,N} (1) 
V N 

works for an arbitrary N, not necessarily prime. A Fourier matrix is unitary, 
but a rescaled matrix Fn = ^/NF' N belongs to H(N, N). In the following we are 
going to use Hadamard matrices with unimodular entries, but for convenience 
we shall also refer to F/v as Fourier matrix. 

For to = 2 it is not difficult [5] to construct matrices H(p,2p). In general, 
the problem of finding all pairs {q, N} for which Butson-type matrices H(q, N) 
do exist, remains unsolved |15|. even though some results on non-existence are 
available |16|. The set of p-th roots of unity forms a finite group and it is 
possible to generalize the notion of Butson-type matrices for any finite group 

DU EU nu. 

In this work we will be interested in a more general case of complex 
Hadamard matrices, for which there are no restrictions on phases of each 
entry of H. Such a matrix, also called biunitary [23], corresponds to taking 
the limit q —> oo in the definition of Butson. In this case there is no analogue 
of Hadamard conjecture, since complex Hadamard matrices exists for any di- 
mension. However, in spite of many years of research, the problem of finding 
all complex Hadamard matrices of a given size N is still open |211 12"2~] . 

Note that this problem may be considered as a particular case of a more 
general issue of specifying all unitary matrices such that their squared moduli 
give a fixed doubly stochastic matrix |23l I24| . This very problem was intensively 
studied by high energy physicists investigating the parity violation and analyzing 
the Cabibbo-Kobayashi-Maskawa matrices |251 1261 1271 |2"S| . 

On one hand, the search for complex Hadamard matrices is closely related to 
various mathcmatical problems, including construction of some *-subalgebras in 
finite von Neumann algebras (2312311312]) analyzing bi-unimodular sequences 
or finding cyclic n-roots [32123 an d equiangular lines [32] • Complex Hadamard 
matrices were used to construct error correcting codes |15| and to investigate 
the spectral sets and Fuglede's conjecture [351 lÏÏül IÏÏ71ISH] . 

On the other hand, Hadamard matrices find numcrous applications in several 
problems of thcoretical physics. For instance, Hadamard matrices (rescaled 

1 These matrices were also called complex Hadamard matrices llülll llll2lll,j| . 
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by 1/y/N to achieve unitarity), are known in quantum òptics as symmetric 
multiports 0Ü] (and are sometimes called Zeilinger matrices) and may be 
used to construct spin models [31] , or schemes of selective coupling of a multi 
qubit system 

Complex Hadamard matrices play a crucial role in the theory of quantum 
information as shown in a seminal paper of Werner 43. They are used in 
solving the Mean King Problem |441 1451 I4ü| . and in finding 'quantum designs' 
47 . Furthcrmore, they allow one to construct 

a) Bases of unitary operators, i.e. the set of mutually orthogonal 
unitary operators, {Uk}^ =1 such that Uk € U(N) and TilllUí = NS^i for 
k,l = l,...,N 2 , 

b) Bases of maximally entangled states, i.e. the set {l^fc)}^! such 
that each \^k) belongs to a composed Hilbert space with the partial trace 
Trjvd^jfeX^jfej) = 1 / N > and the y are mutually orthogonal, (*k|*ü) = 6 k i 0H|, 

c) Unitary depolarisers, i.e. the set {Uk\k=i such that for any bounded 
linear operator A the property Ylk=i ^lAlIk = NTrA 1 holds. 

The problems a)-c) are equivalent in the sense that given a solution to one 
problem one can find a solution to the other one, as well as a corresponding 
scheme of teleportation or dense coding |43j . In particular Hadamard matrices 
are usefull to construct a special class of unitary bases of a group type, also 
called 'nice error basis' |491 15Ü| 

Another application of Hadamard matrices is related to quantum tomog- 
raphy: To determine all N 2 — 1 parameters characterizing a density matrix 
of size N one needs to perform k > N + 1 orthogonal measurements. Each 
measurement can be specified by an orthogonal basis <&„ = {l^^^iLi set for 
u = 1, . . . , k. Precision of such a measurement scheme is optimal if the bases 
are mutually unbiased, i.e. they are such that 

l<4^ s) >l 2 = ^(i-M + <^% . (2) 

If the dimension N is prime or a power of prime the number of maximally 
unbiased bases (MUB) is equal to N + 1 |511 152) . but for other dimensions 
the answer to this question is still unknown [23 02] |SS] . The task of finding 
(k + 1) MUBs is equivalent to finding a collection of k mutually unbiased 
Hadamards (MUH), 

{HieH N }ti- -^=HjH 3 eH N , i>j = l,...,k-l, (3) 

since the set { 1 , H\/ v~N, . . . , Hu /VN} forms a set of MUBs. Here Ti. n denotes 
the set of complex 2 Hadamard matrices of size N. 

The aim of this work is to review properties of complex Hadamard matrices 
and to provide a handy collection of these matrices of size ranging from 2 to 

2 Similarly, knowing unbiased real Hadamard matrices one constructs real MUBs 56 . 
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16. Not only we list concrete Hadamard matrices, the existence of which follows 
from recent papers by Haagerup [21] and Di1;à |22], but also we present several 
other Hadamard matrices which have not appeared in the literature so far. 

2 Equivalent Hadamard matrices and the de- 
phased form 

We shall start this section providing some formal definitions. 

Definition 2.1 A square matrix H of size N consisting of unimodular entries, 
\Hij \ = 1, is called a Hadamard matrix if 

HHÏ = N 1 , (4) 

where < denotes the Hermitian transpose. One distinguishes 

a) real Hadamard matrices, £ R 7 for i,j = 1, . . . , N, 

b) Hadamard matrices of Butson type H(q,N), for which (Hij) q = 1, 

c) complex Hadamard matrices, Hij G C. 

The set of all Hadamard matrices of size N will be denoted by TÍn- 

Definition 2.2 Two Hadamard matrices H\ and H2 are called equivalent, 

written H\ ~ H2 , if there exist diagonal unitary matrices D\ and D2 and per- 
mutations matrices P\ and P2 such that I21f 

Hi = D X P X H 2 P 2 D 2 . (5) 

This equivalence relation may be considered as a generalization of the Hadamard 
equivalence in the set of real Hadamard matrices, in which permutations and 
negations of rows and columns are allowed. 3 

Definition 2.3 A complex Hadamard matrix is called dephased when the en- 
tries of its first row and column are all equal to unity^ , 

Hx,i = H lA = 1 for i = 1, . . . , N . (6) 

Remark 2.4 For any complex N x N Hadamard matrix H there exist uniquely 
determined diagonal unitary matrices, D r = diag(iín, H21, ■ ■ ■ , -ff/vi), and D c = 
diag(l, H\\H\2i ■ ■ ■ ,HhHin), such that [D c ] l 1 — 1 and 

D r ■ H ■ D c (7) 

is dephased. 

3 Such an equivalence relation may be extended to include also transposition and complex 
conjugation 1281 . Since the transposition of a matrix is not realizable in physical systems we 
prefer to stick to the original definition of equivalence. 

4 in case of real Hadamard matrices such a form is called normalised. 
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Two Hadamard matrices with the same dephased form are equivalent. Thus 
the relevant information on a Hadamard matrix is carried by the lower right 
submatrix of size N — 1, called the core [S]. 

It is often useful to define a log Hadamard matrix $, such that 



H k i = e 



i 'i'/. 



(8) 



is Hadamard. The phases entering a log-Hadamard matrix may be chosen 
to belong to [0,27r). This choice of phases implies that the matrix q&/2ir cor- 
responding to a Hadamard matrix of the Butson type H(q, N) consists of zeros 
and integers smaller than q. All the entries of the first row and column of a log 
Hadamard matrix $, corresponding to a dephased Hadamard matrix, are equal 
to zero. 

To illustrate the procedure of dephasing consider the Fourier-like matrix 



of size four, 



Fa 



e ijfe27r/4 w h ere j j. e {1,2,3,4}. Due to this choice 



of entries the matrix F4 is not dephased, but after operation it takes the 
dephased form Fi, 



Fi = 



i -1 

-1 1 

-i -1 

1 1 



-i 1 

-1 1 

i 1 

1 1 



F 4 = 



1 

-1 
1 

-1 



1 

—i 
-1 

i 



(9) 



The corresponding log-Hadamard matrices read 



$4 = 



2tt 



12 3 

2 2 

3 2 10 




$4 = 



2tt 







2 3 
2 
2 1 



(10) 



Note that in this case dephasing is equivalent to certain permutation of rows and 
columns, but in general the role of both operations is different. It is straight- 
forward to perform J7J) which brings any Hadamard matrix into the dephased 
form. However, having two matrices in such a form it might not be easy to 
verify, whether there exist permutations P\ and P2 necessary to establish the 
equivalence relation 101 . This problem becomes difficult for larger matrix sizes, 
since the number of possible permutations grows as NI. What is more, com- 
bined multiplication by unitary diagonal and permutation matrices may still be 
necessary. 

To find necessary conditions for equivalence one may introduce various in- 
variants of operations allowed in (JSJ), and compare them for two matrices in- 
vestigated. In the case of real Hadamard matrices and permutations only, the 
number of negative elements in the dephased form may serve for this purpose. 
Some more advanced methods for detecting incquivalence were recently pro- 
posed by Fang and Ge |S7]. It has been known for many years that for N = 4, 8 
and 12 all real Hadamard matrices are equivalent, while the number of equiva- 
lence classes for N = 16, 20, 24 and 28 is equal to 5, 3, 60 and 487, respectively. 
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For higher dimensions this number grows dramatically: For N — 32 and 36 the 
mimbcr of incquivalent matrices is not smaller than 3,578,006 and 4, 745,357, 
but the problem of enumerating all of them remains open - see e.g. (58| . 

To characterize a complex Hadamard matrix H let us define a set of coefh- 
cients, 

A={A tl :=H ij H^H kl H~r- H = (i, 3, k, l) G {1, . . . ,iV} x4 } , (11) 

where no summation over repeating indeces is assumed and fx stands for a com- 
posite index (i,j,k,l). Due to complex conjugation of the even factors in the 
above definition any concrete value of is invariant with respect to multipli- 
cation by diagonal unitary matrices. 5 Although this value may be altered by a 
permutation, the entire set A is invariant with respect to operations allowed in 
(01. This fact allows us to state the Haagerup condition for inequivalence \'21\ . 

Lemma 2.5 Iftwo Hadamard matrices have different sets A of invariants \TTjj. 
they are not equivalent. 

The above criterion works in one direction only. For example, any Hadamard 
matrix H and its transpose, H T , possess the same sets A, but they need not to 
be permutation equivalent. 

Some other equivalence criteria are dedicated to certain special cases. The 
tensor product of two Hadamard matrices is also a Hadamard matrix: 

H x G H M and H 2 E H N =>> H l ®H 2 <E H M n (12) 

and this fact will be used to construct Hadamard matrices of a larger size. Of 
particular importance are tensor produets of Fourier matrices, Fm®Fn £ TLmn ■ 
For arbitrary dimensions both tensor produets are equivalent, Fn®Fm — Fm® 
Fj\f. However, their equivalence with Fmn depends on the number theoretic 
property of the product M ■ N : the equivalence holds if M and N are relatively 
prime 

To classify and compare various Hadamard matrices we are going to use 
the dephased form (JüJ) thus fixing the first row and column in each compared 
matrix. However, the freedom of permutation inside the remaining submatrix 
of size (N — 1) does not allow us to specify a unique 'canonical form' for a 
given complex Hadamard matrix. In other words we are going to list only 
certain representatives of each equivalence class known, but the reader has to 
be warned that other choices of equivalent representatives arc equally legitimate. 
For instance, the one parameter family of N = 4 complex Hadamard matrices 
presented in contains all parameter-dependent elements of H in its lower- 

right còrner, while our choicc Ij65(l with variable phases in second and fourth 
row is due to the fact that such an orbit stems from the Fourier matrix F^. 

Identifying matrices equivalent with respect to JHJ we denote by 
Gn = U N / ~ (13) 

5 Such invariants were used by physicists investigating unitary Kobayashi-Maskawa matri- 
ces of size 3 and 4 |2S]H3I2ZII2ID· 
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every set of representatives of different equivalence classes. 

Interestingly Qn is known only for N — 2, 3, 4, 5, while the problem of finding 
all complex Hadamard matrices for N > 6 remains unsolved. In particular, 
compiling our list of Hadamard matrices we took into account all continuous 
families known to us, but in several cases it is not clear, whether there exist any 
equivalence relations between them. 

3 Isolated Hadamard matrices and continuous 
orbits of inequi valent matrices 

In this and the following sections we shall use the symbol o to denote the 
Hadamard product of two matrices, 

[H! o H 2 ] t j = [Hi]^ ■ [H 2 ] itj , (14) 

and the EXP symbol to denote the entrywise exp operation on a matrix, 

[EXP(fl)]. i = exp ([R] id ) . (15) 
3.1 Isolated Hadamard matrices 

Definition 3.1 A dephased N x N complex Hadamard matrix H is called iso- 
lated if there is a neighbourhood W around H such that there are no other 
dephased complex Hadamard matrices in W. 

To have a tool useful in determining whether a given dephased complex 
Hadamard matrix of size N is isolated, we introduce the notion of defect: 

Definition 3.2 The defect d{H) of an N x N complex Hadamard matrix H 
is the dimension of the solution space of the real linear system with respect to a 
matrix variable R € ïïí N 

Rij =0 j€{2,...,N} (16) 

Ri,i =0 ie{i,...,N} (17) 

N 

J2 H i,kH j , k (Ri,k-R j ,k) - l<i<j<N (18) 

k—1 

The defect allows us to formulate a one way criterion: 

Lemma 3.3 A dephased complex Hadamard matrix H is isolated if the defect 
of H is equal to zero. 

If N is prime then the defect of Fm is zero, so the Fourier matrix is isolated, 
as earlier shown in |20l I62| . For any composed N the defect of the Fourier 
matrix is positive. For instance, if the dimension N is a product of two distinct 
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primes, then d(F pq ) = 2(p — l)(q — 1), while for powers of a prime, N — p k with 
k > 2, the defect reads 

d(F pk ) = p k - 1 [k(p-l)-p} + l . (19) 

An explícit formula for d(F^r) for an arbitrary composed N is derived elsewhere 
63 . In that case the Fourier matrix belongs to a continuous family, as examplcs 
show in scction[S] 

The reasoning behind the criterion of Lemma 13.31 runs as follows: 

Any dephased complex Hadamard matrices, in particular those in a neigh- 
bourhood of a dephased complex Hadamard matrix H, must be of the form: 

HoEXP{i-R) (20) 

where an n x n real matrix R satisfies 'dephased property' and unitarity condi- 
tions for lt2U|l : 



JY 



Rí,j 


= 


3 e {2, . . 


.,N} 


(21) 


Ri,i 


= 


i e {i,... 


,N} 


(22) 


-Rj.k) 


= 


í<i<j 


< N 


(23) 



k=l 

We will rewrite these conditions using a real vector function /, whose co- 
ordinate functions will be indexed by the vàlues (symbol sequences) from the 
set X, related to the Standard index set by a fixed bijection (one to one map) 
13: 1 — > {1, 2, . . . , (2N - 1) + (A^ 2 - N)}. The set 1 reads: 

1 = { (1,2), (1,3), ... , (1,N) } U { (1,1), (2,1), ... , (A,l) } 

U { (i,j,t); l<i<j<N and í e {Re, Im} } (24) 

For simplicity of notation we further write fi, i € X to denote fp(i)- 

Similarly, a fixed bijection a : {1, . . . , N} x {1, . . . , N} — > {1, . . . , N 2 } 
allows matrix indexing the components of a real N 2 element vector variable R, 
an argument to /, and we write Rk,i to denote R a (k,i) 

The (2N — 1) + (N 2 — N) element function vector / is defined by the formulas: 

/(ij) = R i,i 

f(i,l) = RiS 

/ A \ 

i(fl,. fe -fl 3 , fc ) 



/(.j,im) = Im -i ■ H hk H j:k e l{ 



J G {2, . . 


.,N} 


(25) 


ÍG{1,... 


,N} 


(26) 


l<i<j 


< N 


(27) 


1 < i < 3 


< N 


(28) 



fe=l 
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Conditions I|21I22I23I) can now be rewritten as: 



f(R) = (29) 

where the a(i,j)-th coordinate of a real variable vector R represents the i, j-th 
entry of the corresponding matrix R sitting in (|20|l . 

The value of the linear map Df : R" 2 — ► R (2N-i)+(A^i) ^ being thc 
differential of / at 0, at R is the vector: 

[Df (R)] (1J) = R hj je{2,...,N} (30) 
[Df (R)} {il) = Ri A ie{l,...,N} (31) 



l<i<j<N (32) 



P/oí^lftj» = Im \J2 H i,kHj,k (Ri,k ~ R j>k )\ l<i<j<N (33) 

where again indexing for / defined by (3 is used. 

It is clear now that the kernel of the differential, {R G E* 2 : Df (R) = 0}, 
corresponds to the solution space of system (|16I17I18|I . in which R now takes 
the meaning of an input variable vector to /, with indexing determined by a 

Note that the N 2 — N equation subsystem (JTHJ: 

N 

J2 H i,kH j ,k(Ri,k-Rj,k) = 1 < i < i < iV (34) 

k=i 

is solved at least by the (2N — l)-dimensional real space spanned by 2N — 1 
vectors, defined by: 

Mi f0t(M) o«^; W1XÜ1 ' ^ÍV.., W} (35, 

and 

fi fei = í i for(M)e«x{i,...,Ar} 

\ otherwise ' 1 ' ' J v ; 

that is, by vectors, if treated as matrices (with the i, j-th entry being equal to 
the a(i, j)-th coordinate of the corresponding variable vector), forming matrices 
with either a row or a column fillcd all with l's and the other entries being O's. 

If the defect of H equals 0, then the overall system (|16ll7ll8fl is solved only 
by 0, the differential Df has full rank N 2 , i.e. dim ^_D/ (R Aí2 )^) = iV 2 , and we 
can choose an N 2 equation subsystem: 

f(R) = (37) 



9 



of 1(23 such that thc diffcrcntial Df at is of rank N 2 , i.e. dim \Df (R N2 )J = 

N 2 , and thus / satisfies The Inverse Function Theorem. The theorem implics 
in our case that, in a neighbourhood of 0, R = is the only solution to (|37() . 
as well as the only solution to 1(25} . 

Also in this case, let us consider the differential, at 0, of the partial function 
vector / , given by 1(27128} . This differential value at R, Df u (R), is given 
by the partial vector I|32I33[) . Since the 'remaining' differential, corresponding 
to the dephased 'property condition' part of /, defined by (|30I31|I . is of rank 
(2JV-1), the rank of Df u is equal to N 2 -(2N-l). Recali that from considering 
above the minimal solution space of system <|18ll . it cannot be greater than 
TV 2 — (2N — 1). Were it smaller, the rank of Df would be smaller than N 2 , 
which cannot be if the defect of H is 0, see above. 

Then one can choose an N 2 — (2N — 1) equation subsystem f u (R) = 
of system f u {R) — 0, with the full rank: 

dim (Df u (ïïl N2 )^ = N 2 — (2N - 1) (38) 

thus defining a (2N — 1) dimensional manifold around 0. This manifold gener- 
ates, by Q20[l. the (2N — 1) dimensional manifold containing all, not necessarily 
dephased, complex Hadamard matrices in a neighbourhood of H. In fact, the 
latter manifold is equal, around H, to the (2N — 1) dimensional manifold of 
matrices obtained by left and right multiplication of H by unitary diagonal 
matrices: 

{ diag(e iai ,...,e ÍQjv ) H ■ diag(l, e i& , . . . , e i/3jv ) } (39) 

3.2 Continuous orbits of Hadamard matrices 

The set of inequivalent Hadamard matrices is finite for N = 2 and N = 3, but 
already for N = 4 there exists a continuous, one parameter family of equivalence 
classes. To characterize such orbits we will introduce the notion of an afHne 
Hadamard family. 

Definition 3.4 An affine Hadamard family H(7Z) stemming from a de- 
phased N x N complex Hadamard matrix H is the set of matrices satisfying 
, associated with a subspace 1Z of a space of all real N x N matrices with 
zeros in the first row and column, 

H(n) = {Ho EXP(i -R) : R ell} . (40) 



The words 'family' and 'orbit' denote submanifolds of M. 2n consisting purely 
of dephased complex Hadamard matrices. We will often write H(a,\, . . . , a m ) if 
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1Z is known to be an m-dimensional space with basis R\, . . . , R rn . In this case, 
H(a%, . . . , a m ) will also denote the element of an affinc Hadamard family: 



H(a u ..., 



H(R) = H o EXP(i • R) 



(41) 



where R = a\ ■ R\ + ... + a, 



Ri 



rn 



An affinc Hadamard family H(1Z) stemming from a dephased NxN complex 
Hadamard matrix H is called a maximal affine Hadamard family when it 
is not contained in any larger affine Hadamard family H(1Z') stemming from H: 



Calculation of the defect of H, defined in the previous section, is a step 
towards determination of affine Hadamard families stemming from H: 

Lemma 3.5 There are no affine Hadamard famílies stemming from a dephased 
NxN complex Hadamard matrix H if it is isolated, in particular if the defect 
of H is equal to 0. 

Lemma 3.6 The dimension of a continuous Hadamard orbit stemming from a 
dephased Hadamard matrix H is not greater than the defect, d(H) . 

A lower bound d c (N) for the maximal dimensionality of a continuous orbit 
of inequivalent Hadamard matrices of size N was derived by Ditjà |22|. Inter- 
estingly, for powers of a prime, N = p k , this bound coincides with the defect 
(|19f) calculated at the Fourier matrix, which provides an upper bound for the 
dimension of an orbit stemming from F^. Thus in this very case the problem 
of determining d c (N) restricted to orbits including Fn is solved and we know 
that the maximal affine Hadamard family stemming from Fn is not contained 
in any non-affine orbit of a larger dimension. 

Finally, we introduce two notions of equivalence between affine Hadamard 
families: 

Definition 3.7 Two affine Hadamard families stemming from dephased NxN 
complex Hadamard matrices H\ and H^: H±(TZ') and Hi^R."), associated with 
real matrix spaces 1Z' and 1Z" of the same dimension, are called permutation 
equivalent if 

there exist two permutation matrices P r and P c such that 



i. e. there is one to one correspondence, by row and column permutation, between 
the elements of Hi(lZ') and H2{1Z"). 




(42) 



H 2 {n") = p r -H x {n')-p c 



(43) 
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Note that permutation matrices considered in the above definition must not 
shift thc first row or column of a matrix. 

Definition 3.8 Two affine Hadamard famílies H±(Tt') and H^iTZ!') stemming 
from dephased N x N complex Hadamard matrices H\ and H2, associated with 
real matrix spaces 1Z' and 1Z" of the same dimension, are called cognate if 

VBeH 2 (K") BAeH^TZ') B ~ A T , (44) 
VA € ffi(Te') 3B e H 2 (n") A ~ B T . (45) 

The family H(1Z) is called self cognate if 

VBeH(K) 3AeH(K) B~A t (46) 



4 Construction of Hadamard matrices 

4.1 Same matrix size: reordering of entries and conjuga- 
tion 

If H is a dephased Hadamard matrix, so are its transpose H T , the conjugated 
matrix H and Hermitian transpose . It is not at all obvious, whether any 
pair of these is an equivalent pair. However, in some special cases it is so. 

For example, the dephased forms of a Hadamard circulant matrix C and its 
transpose C T are equivalent since C and C T are always permutation equivalent 
(see also the remark on the top of p.319 in [21]): 

C T = P T CP (47) 

where C is an N x iV circulant matrix Cjj = a^-j m od N f° r some x £ (D , and 
P = [ei, eAr, ejv_i, . . . ,e^[, where are the Standard basis column vectors. 

On the other hand, there are infinitely many examples of equivalent and 
inequivalent pairs of circulant Hadamard matrices C, C, the same applying to 
their dephased forms. 

Apart from transposition and conjugation, for certain dimensions there exist 
other matrix reorderings that preserve the Hadamard structure. Such operations 
that switch substructures of real Hadamard matrices to generate inequivalent 
matrices have recently been discussed in |SHj. It is likely that these methods 
may be useful to get inequivalent complex Hadamard matrices. In this way only 
a finite number of inequivalent matrices can be obtained. 

4.2 Same matrix size: linear variation of phases 

Starting from a given Hadamard matrix H in the dephased form one may inves- 
tigate, whether it is possible to perform infinitesimal changes of some of (N— l) 2 
phases of the core of H to preserve unitarity. Assuming that all these phases 
(<í>fe;, from Eq. (JHJ fe, l = 2, . . . , N) vary linearly with free parameters one can 
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find analytical form of such orbits, i.e. affine Hadamard famílies, stemming from 
e.g. Fourier matrices of composite dimensions |60| . 

To obtain affine Hadamard families stemming from H, one has to consider 
all pairs of rows of H. Now, taking the inner product of the rows in the i, j-th 
pair (1 < i < j < N), one gets zero as the sum of the corresponding vàlues in 
the sequence: 

( //,.; -77,.: , II, ,> ■ 77,., , ... . H iíN -Hj, N ) (48) 

Such sequences will further be callcd chains, thcir subsequences - subchains. 
Thus l|48l) features the i, j-th chain of H. A chain (subchain) is closed if its 
elements add up to zero. As (l/VTV) • H is unitary, all its chains are closed. 
It is not obvious, however, that any of these chains contain closed subchains. 

Let us now construct a closed subchain pattern for H . For each chain of 
H, let us split it, disjointly, into closed subchains: 

H. m ■ H. , (,) , H , (a) • H (s) , . . . , H ( s) • H ( ») 

L -!^\ J-^i 6 )"'2 Jf"·2 ' ( s ) (s) 



(49) 



where 



s 6 {1, 2, . . . , rjj} designates subchains, r^j being the number of closed 
subchains the i, j-th chain is split into 

Ul=i{ fc í S) > ■ ■ ■ , k p Í)} = {1, 2, . . . , N}, being the length of the s-th sub- 

chain 

{fc^\...,^}n{fc 1 - ) ,...,^ 2 ) ) } = if Sl ^.s 2 

and splitting of {1, 2, . . . , N} into {k[ , ■ ■ ■ , fe^j } is done independently for each 
chain of H, that is the /c-values above in fact depend also on i,j. 

A pattern according to which all the chains of H are split, in the above way, 
into closed subchains, will be called a closed subchain pattern. 

A closed subchain pattern may give rise to the affine Hadamard family, 
stemming from H , corresponding to this pattern. The space 1Z, generating this 
H(7Z) family (see Definition l3.4|) . is defined by the equations: 

R ltj = je{2,...,N}, (50) 
Ri,l = t G {!,..., N}, (51) 

(as H(1Z) is made of dephased Hadamard matrices) and, for all 1 < i < j < N 
and s £ {1, 2, . . . , r iyj } , 



= = IM>-M)' ' (52) 
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where the sets of indices {k[ s , . . . , k (j-, } correspond to the considered pattern 
splitting of the i,j-th chain into closed subchains l|49|l . Recali that the sets of 
fc's depend on i, j's, ommited for simplicity of notation. 

If the system (150151152(1 yields a nonzero space TZ, then H(TZ) is the affine 
Hadamard family corresponding to the chosen closed subchain pattern. Actu- 
ally, any affine Hadamard family stemming from H is generated by some space 
TZ contained in a (probably larger) space TZ', corresponding to some closed sub- 
chain pattern for H. The respective theorem and its proof will be published in 

It may happen that the system 1(51)151152(1 . for pattern Pi^, defines space 
TZ, which is also obtained as the solution to (|50I51I52() system shaped by 
another pattern P2 < " H \ imposing stronger conditions of type 1)52(1. as a result 
of there being longer subchains in P2^ íf ·' composed of more than one subchain 
of P].^, for a given pair i,j. If this is not the case, we say that TZ is strictly 
associated with pattern P±^ H \ Then the subchains of the i,j-th chain of 
H(R) = Ho EXP(i ■ R), ReTZ, distinguished according to pattern Pi (H) : 

^i.itS'^i.fcW exp (* ■ ( R i,k^ - i? J ,fe< s) )) ' 
H („ H w exp ( i • ( R w - R w )) ) (53) 

rotate independently as i? runs along 7?. satisfying 1(52(1 . 

Maximal affine Hadamard families stemming from H are generated by maxi- 
mal, in the sense of l(42|l , solutions to I(5UI51I52|) systems shaped by some specific 
closed subchain patterns, which we also call maximal. We have been able to 
find all of these for almost every complex Hadamard matrix considered in our 
catalogue. To our undcrstanding, however, it becomes a serious combinatorial 
problem already for N — 12. For example, for a real 12 x 12 Hadamard matrix, 
each chain can be split in (6!) 2 ways into two element closed subchains only. 

Fortunately, as far as Fourier matrices are concerned, the allowed maximal 
subchain patterns are especially regular. We are thus able to obtain all maximal 
affine Hadamard families stemming from Fjy for an arbitrary N (6U( , as we have 
done for TV < 16. 

An alternative method of constructing affine Hadamard families, developed 
by DfÇà [22], is presented in section |4~51 We also refer the reader to the article 
by Nicoara |2Ü( , in which conditions are given for the cxistcncc of one parame- 
ter families of commuting squares of finite dimensional von Neumann algebras. 
These conditions can be used to establish the existence of one parameter fam- 
ilies of complex Hadamard matrices, stemming from some H, which are not 
assumed to be affine. 
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4.3 Duplication of the matrix size 



Certain ways of construction of real Hadamard matrices of an extended size 
work also in the complex case. If A and B belong to TÍn then 



H 



A B 

A -B 



€ TÍ2N 



(54) 



Furthermore, if A and B are taken to be in the dephased form, so is H . 
This method, originally due to Hadamard, can be generalized by realizing 
that B can be multiplied at left by an arbitrary diagonal unitary matrix, 
E = diag(l, e ldl , . . . , e 1 '*™ -1 ). If A and B depend on a and b free parameters, 
respectively, then 

A EB ~ 
A -EB 



H' = 



(55) 



represents an (a + b + N - 
2N in the dephased form. 



l)-parameter family of Hadamard matrices of size 



4.4 Quadruplication of the matrix size 

In analogy to Eq. I|54|l one may quàdruple the matrix size, in a construction 
similar to that derived by Williamson j^I] from quaternions. If A, B,C, D G Hn 
are in the dephased form then 

' A B C D 

A -B C -D 
A B -C -D 
A -B -C D 

and it has the dephased form. This form is preserved, if the blocks B, C and 
D are multiplied by diagonal unitary matrices E\ , E2 and E3 respectively, each 
containing unity and N — 1 free phases. Therefore l|56[) describes an [a + b + 
c + d+'ò(N — l)]-dimensional family of Hadamard matrices [22], where a,b,c,d 
denote the number of free parameters contained in A, B, C, D, respectively. 



€ TÍ4N 



(56) 



4.5 Generalized method related to tensor product 

It is not difficult to design a similar method which increases the size of a 
Hadamard matrix by eight, but more generally, we can use tensor product 6 
to increase the size of the matrix K times. For any two Hadamard matrices, 
A € Hk and B S TÍm, their tensor product A <£> B S TLkm- A more general 
construction by Dit,à allows to use entire set of K (possibly different) Hadamard 
matrices {-Bi, . . . , Bk } of size M. Then the matrix 





' A 1X B X 


A12E2B2 


AirEkBk 


H = 










AkxB x 


AK2E2B2 ■ 


■ AkkEkBk 



Tensor products of Hadamard matrices of the Butson type were investigated in |39|. 
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of size N — KM is Hadamard |22|. As before we introduce additional free 
phases by using K — 1 diagonal unitary matrices Ek ■ Each matrix depends on 
M — 1 phases, since the constraint [Ek] 1 1 = 1 for k = 2, . . . , K is necessary to 
preserve the dephased form of H . Thus this orbit of (mostly) not equivalent 
Hadamard matrices depends on 

K 

d = a + J2 b ] + (K-1)(M-Ï) (58) 
i=i 

free parameters. Here a denotes the number of free parameters in A, while bj 
denotes the number of free parameters in Bj. A similar construction giving 
at least (K — 1)(M — 1) free parameters was given by Haagerup I n the 

simplest case K ■ M = 2 • 2 these methods give the Standard 1-parameter N = 4 
family i|t)5|) . while for K ■ M = 2-3 one arrives with (2 — 1)(3 — 1) = 2 parameter 
family (|69() of (mostly) inequi valent N = 6 Hadamard matrices. 

The tensor product construction can work only for composite N, so it was 
conjectured |5Ü] that for a prime dimension TV there exist only fmitely many 
inequivalent complex Hadamard matrices. However, this occurred to be false 
after a discovery by Petrescu |62| . who found continuous families of complex 
Hadamard matrices for certain prime dimensions. We are going to present his 
solution for N — 7 and N = 13, while a similar construction jïï2] works also 
for N — 19, 31 and 79. For all primes N > 7 there exist at least three isolatcd 
complex Hadamard matrices, see [21], Thm. 3.10, p.320. 

5 Catalogue of complex Hadamard matrices 

In this section we list complex Hadamard matrices known to us. To save space 
we do not describe their construction but present a short characterization of 
each case. 

Each entry H , provided in a dephased form, represents a continuous (2A^ — 
l)-dimensional family of matrices H obtained by multiplication by diagonal 
unitary matrices: 

H(ai, ...,oín,02, —,Pn) = Di{a u . . .,a N ) ■ H ■ D 2 {f3 2 , ■ ■ ■ ,0n) (59) 

where 

D 1 (ai,... ) a J v) = diag(e ia \...,e ÍQ ") 

D 2 (Jh,...,0N) = diag(l, e ife ,...,e i/3N ) 

Furthermore, each H in the list represents matrices obtained by discrete 
permutations of its rows and columns, and equivalent in the sense of J5J- For 
a given size N we enumerate families by capital letters associated with a given 
construction. The superscript in brackets denotes the dimension of an orbit. For 

(2) 

instance, Fq represents the two-paramctcr family of N — 6 complex Hadamard 
matrices stemming from the Fourier matrix Fq. Displaying continuous families 
of Hadamard matrices we shall use the symbol • to denote zeros in phase vari- 
ation matrices. For completeness we shall start with a trivial case. 
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5.1 JV = 1 

f[ 0) =F 1 = [1]. 



5.2 N = 2 

All complex 2x2 Hadamard matrices are equivalent to the Fourier matrix F% 
(|2T], Prop.2.1, p.298): 



F 2 =F 2 



(o) 



1 1 

1 -1 



(60) 



The complex Hadamard matrix F% is isolated. 



The set of inequivalent Hadamard matrices of size 2 contains one element 
only, Q 2 = {F 2 }. 



5.3 N = 3 

All complex 3x3 Hadamard matrices are equivalent to the Fourier matrix F% 
(j2U, Prop.2.1, p.298): 



F 3 = = 



1 1 



1 



1 w w 
1 w w 



(61) 



where w = exp(i • 27r/3), so w 3 = 1. 

i*3 is an isolated complex Hadamard matrix. 



The set of dephased representatives can be taken as Ç/3 = {-F3}. 



5.4 N = A 



Every 4x4 complex Hadamard matrix is equivalent to a matrix belonging to 
the only maximal afiinc Hadamard family Í4 (<z) stemming from F4 |21| . The 



F^ (a) family is given by the formula: 



F 4 (1) (a) =F 4 oEXP(i-iï p(1) (a)) 



(62) 



where 



Fi 



1 


1 


1 


1 




" 1 


1 


1 


1 


1 


w 


2 

w 


w 3 




1 


i 


-1 


i 


1 




1 


2 




1 


-1 


1 


-1 


1 


w 3 




w 






1 


i 


-1 


i 



(63) 
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where w — cxp(i • 27r/4) = i, so u> 4 = 


1, w 2 




-1, and 






• • 


• 


• 






(a) = 


• a 


• 


a 




• • 

• a 


• 
• 


• 

a 




Thus the orbit stcmming from F4 reads: 










" 1 


1 


1 




1 




1 i 

1 


• e ia 

-1 


-1 
1 


— i 


• e ia 
-1 




1 - 


i-e ia 


-1 


i ■ 





(64) 



(65) 



The above orbit is permutation equivalent to: 





F 2 


[ía]i,2 


D(a) 


F 2 ' 




F 2 


[F2} 2 , 2 


D(a) 


F 2 \ 



(66) 



where D(a) is the 2x2 diagonal matrix diag(l, e 1Q ). 

This orbit is constructed with the Dinà's method [22J , by setting K = M 
2, A = B X =B 2 =F 2 , E 2 =D{a) in (52|. 



It passes through a permuted Í4: 

F 4 = 5 4 (1) (tt/2) • [ei,e3,e2,e 4 ] T 



(67) 



where is the i-th Standard basis column vector, and through F 2 ®F 2 = F^ (0). 
Note that F4 and F 2 ® F2 are, according to [Sü|, inequivalent. 

The Fj orbit of (|62l) is symmetric, so it is self-cognate. Replacing a by 
a + 7r yields F4 (a) with the 2-nd and 4-th column exchanged, so F4 (a) ~ 
F^ia + Tr). 

The set of dephased representatives can be taken as Gi — {Fj (<z) : a £ 

[0.7T)} 



5.5 iV = 5 

As shown by Haagerup in |23 (Th. 2.2, p. 298) for N = 5 all complex Hadamard 
matrices are equivalent to the Fourier matrix F5: 



F 5 = Ff = 



1 


1 


1 


1 


1 


1 


w 




w z 




4 

w 


1 


w 2 


«> 4 


U) 


w 3 


1 


w 3 


w 


4 

tir 


2 
ur 


1 


4 

w 


w 3 


ur 


w 



(68) 
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where w = exp(i • 27r/5) so w 5 = 1. 
The above matrix i*g is isolated. 



The set of dephased representatives can be taken as G5 = {F5}. 



5.6 N = 6 

5.6.1 Orbits stemming from F§ 

The only maximal affine Hadamard families stemming from Fq are: 
F^\a,b) = F 6 oEXP (i- R F w(a,b) 

T 



= F R o EXP i 



R F m(a,b) 

1 ft 



where 





' 1 


1 


1 


1 


1 


1 




1 


w 


2 
tir 




4 

UT 


w 5 




1 


2 

w 


4 

w 


1 


2 
tíT 


4 

w 


1 




1 




1 


w 3 




1 


4 

i/> 


2 

w 


1 


4 
UT 


2 

w 




1 






w 3 




w 


2tt/6) so 


w 6 = 


= 1, 


w 3 = 


— 1, and 










• a 


6 


• a 


b 


-R p (2) ( 


0,6) 




• a 


ò 


• a 


b 


• a 


6 


• a 


b 



(69) 
(70) 



(71) 



(72) 



The above families (|69l) and H70JI are cognate. The defect (see Def. ) reads 



d(F 6 ) = 4. 

At least one of the above orbits is permutation equivalent to one of the orbits 
obtaincd using the Ditjà's mcthod, cithcr 





F 2 


m h2 


D(ai) 


F 2 


[F3 


1,3 


D(a 2 ) 


■F 2 ■ 




F 2 


m 2 , 2 


D(ai) 


F 2 


[F 3 


2,3 


D(a 2 ) 


■F 2 


. i^kx 


F 2 


m,, 2 


D(ai) 


F 2 


[F 3 


3,3 


D(a 2 ) 


■F 2 _ 



(73) 



where D(a) is the 2x2 diagonal matrix diag(l, e 1Q ), by setting K = 3 , M 
2, A = F 3 , B 1 = B 2 = B 3 =F 2 , E 2 = D(oi), E 3 = D(a 2 ) in 
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or 



F^(f3 í ,f3 2 ) = 





[F 2 }!, 2 


D{p! 


P2) 


■F 3 ■ 


■ F 3 


[F2] 2 , 2 


D(Jh 


h) 


■F 3 



(74) 



where D(f3i, (3 2 ) is the 3x3 diagonal matrix diag(l, e 1 ^ 1 , e 1 ^ 2 ), by setting K = 
2, M = 3, A = F a , B 1 = B 2 = F 3 , E 2 = D((3 u f3 2 ) in jE). 

The above statement is true because the orbits (|73[) and l|74|) pass through 



F 6 (2) (0,0) = F 3 ® F 2 and F 6 ^(0,0) = F 2 <g> F 3 , both of which are, accord- 
ing to [Sü], permutation equivalent to Fq. Thus l|73|) and l)74|) are maximal 
afhne Hadamard families stemming from permuted .Fe's. These orbits were 
constructed in the work of Haagerup |21j . 

(2) 

A change of the phase a by tt in the family Fg corresponds to the exchange 

(2) 

of the 2-nd and 5-th column of Fg (a,b), while a change of 6 by tt is equiv- 

(2) 

alent to the exchange of the 3-rd and 6-th column of F 6 (a, b). This implies 
(permutation) equivalence relation: 



M2) 



F 6 (2) (a, 6) ~ F 6 (2) (a + tt, 6) ~ F 6 (2) (a, 6 + tt) ~ F 6 (2) (a + tt, 6 + tt) 



(75) 



5.6.2 1-parameter orbits 

There are precisely five permutation equivalent 1-parameter maximal afHne 
Hadamard families stemming from the symmetric matrix D 6 : 



D fí 



1 

-1 
i 

—i 
—i 
i 



1 

i 

-1 



1 1 

—i i 
i —i 



i -1 i 
-i i -1 
—i —i i 



1 

i 

—i 
—i 
i 

-1 



(76) 



namely 



where 



R (d (c) = 



D 6 o EXP ( i • R D m 



(c) 



• ce» 

— c • • — c 

— c • • — c 

• ce» 



(77) 



(78) 



20 



and 

P 1 ·D 6 1) (c)·P 1 T = D a oEXP(i·P 1 R D M{c)lf y j (79) 

P 2 -D { 6 1) (c)-P? = D 6 oEXP (i -P^a^Pj) (80) 

P 3 -D^ ) (c)-P 3 T = D 6 o-EXP(i-P 3 R DW (c)P£) (81) 



where 



P 4 ·4 1) (c)·P 4 T = D 6 oEXP (i-P4Ü n (i)(c)Pf ) (82) 



Pi = [ei,e3,e 2 ,e 6 ,e 5 ,e4] (83) 

P 2 = [ei,e4,e 3 ,e 2 ,e 6 ,e 5 ] (84) 

P3 = [ei,e 6 ,e 2 ,e 3 ,e4,e 5 ] (85) 

P4 = [ei,e 5 ,e4,e 3 ,e 2 ,e 6 ] (86) 



and e 2 ; denotes the ï-th Standard basis column vector. 
We have the permutation equivalence 



d£\c + tt)=P t -d£\c)-P_ for P_ = [e 1; e 2 , e 3 , e B , e 4 , e 6 ] (87) 



; v<- T /i; — -1 _ • 

Also 



D<£\-c)=(D<£\cj) (88) 
thus {P^ 1} (c) : c G [0,2tt)} and { (c)) T : c € [0,2tt)} are equal sets, that 



is Pg 1 " 1 is a self-cognate family 



None of the matrices of and ÍP@ 2 M of l|69|) and l|7U|l are equivalent to 
" of the matrices of . 

Obviously, the above remarks apply to the remaining orbits stemming from 



D, 



6- 



The D { 6 1} orbit was presented in |22| in the Introduction, and its 'starting 
point' matrix P 6 01 lT76"|) even in the earlier work [51] p.307 (not dephased). 

5.6.3 The 'cyclic 6 roots' matrix 

There exists anothcr, inequi valent to any of the above 6x6 matrices, complex 
Hadamard matrix derived in [21] from the results of 32 on so called cyclic 
6-roots. 



The matrix Cq ^ below is circulant, i.e. it has the structure 



°6 



xçi-j m od 6)+i, where 



x=[l, i/d, -1/d, -i, -d, id] T (89) 
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*= — + '{-) 

is a root of the equation d 2 — (1 — \/3)d +1 = 0. 

The matrix Cq ^ and its dephased form Cg ^ read: 



1 


i d 


-d 


—i 


-d- 1 


irJ- 1 


id- 1 


1 


i d 


-d 


—i 


-d- 1 


-d- 1 


id- 1 


1 


i d 


-d 


—i 


—i 


-d- 1 


id- 1 


1 


i d 


-d 


-d 


—i 


-d- 1 


id- 1 


1 


i d 


i d 


-d 


—i 


-d- 1 


y- 1 


1 



(91) 



1 


1 


1 


1 


1 


1 


1 


-1 


-d 


-d 2 


d 2 


d 


1 


-d- 1 


1 


d 2 


-d 3 


d 2 


1 


-d- 2 


d- 2 


-1 


d 2 


-d 2 


1 


d- 2 


-d" 3 


d- 2 


1 


-d 


1 


d- 1 


d- 2 


-d- 2 


-d- 1 


-1 



(92) 



The circulant structure of Cg implies that it is equivalent to yC 6 j , so 

Cf ^ (Cf) 7 (see 63). Thus C<°> ~ cf ( cf ~ Ó®). 

No affine Hadamard family stems from Cg . However, we do not know if 
the cyclic 6-roots matrix is isolated since the defect d(C^ ] ) = 4 and we cannot 
exclude existence of some other orbit. 



5.6.4 The 'spectral set' 6x6 matrix 

Another complex Hadamard matrix found by Tao (35| plays an important role 
in investigation of spectral sets and disproving the Fuglede's conjecture 
It is a symmetric matrix S'g '', which belongs to the Butson class -ff(3, 6), so its 
entries depend on the third root of unity ui = exp(i • 2w/3), 



1 


1 


1 


1 


1 


1 


1 


1 


UJ 


UJ 


ÜJ 2 


UJ 2 


1 


ÜJ 


1 


ÜJ 2 


UJ 2 


UJ 


1 


ÜJ 


UJ 2 


1 


ÜJ 


UJ 2 


1 


ÜJ 2 


UJ 2 


UJ 


1 


íú 


1 


UJ 2 


U> 


uJ 2 


íú 


1 



(93) 
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Thus the corresponding log-Hadamard matrix reads 

" 

112 2 

1 2 2 1 (94) 
12 12 ■ [y > 

2 2 10 1 
2 12 10 

Its defect is equal to zero, d{S§) = 0, hence the matrix (|9"3"|l is isolated. 
Spectral sets allow to construct certain Hadamard matrices for other composite 
dimensions (see prop. 2.2. in |37]), but it is not yet established, in which cases 
this method yields new solutions. 



$s 6 



2tt 

T 



5.7 N = 7 

5.7.1 Orbits stemming from Fj 

F? is an isolated 7x7 complex Hadamard matrix: 



F 7 = F. 



(0) 



1 


1 


1 


1 


1 


1 


1 


1 


w 


2 

w 


w 3 


4 

W 


w 5 


íc (i 


1 


2 

w 


4 

W 


w 6 


w 


w 3 




1 


w 3 


w 6 


w 


w 5 


w 


4 

W 


1 


4 

w 


w 


w 5 


o 

w 


w 6 


w 3 


1 


w 5 


w 3 


w 


w 6 


4 

W 


w 2 


1 


w 6 


w 5 


4 

W 


w 3 


2 

vr 


w 



(95) 



where w = exp(i • 2w/7), so w 7 = 1. 
5.7.2 1-parameter orbits 

There are precisely three permutation equivalent 1-parameter maximal affinc 
Hadamard families stemming from the symmetric matrix being a permuted 
'starting point' for the 1-parameter or bit found by Petrescu |62|: 



1 


1 


1 


1 


1 


1 


1 


1 


w 


4 

w 


w 5 


w 3 


w 3 


u< 


1 


4 

W 


w 


w 3 


w 5 


w 3 


u< 


1 


w 5 


w 3 


w 


4 

w 


w 


w 3 


1 


w 3 


w 5 


4 

w 


w 


w 


w 3 


1 


w 3 


w 3 


w 


w 


4 
líT 


w 5 


1 


w 


w 


w 3 


w 3 


W 5 


4 

W 


so 


w 6 = 


1, w 


3 


-1. 







(96) 
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They are: 



where 



and 



P 7 (1) (o) = P70 EXP ( i • R (i, (o) 



fipd)(a) 



• a 


a 


• 


• • • 


• a 


a 


• 


• • • 


• • 


• 


— a 


— a • • 


• • 


• 


— a 


—a • • 



P 2 -P$ 1 \a)-P? = P 7 o-EXP(i-P 2 R p w(a)P[ 
where 

Pi = [ei, e 7 , e 3 , e 4 , e 6 , e 5 , e 2 ] P 2 = [ei, e 2 , e 7 , e 6 , e 5 , e 4 , e 3 ] (101) 
and ei denotes the z-th Standard basis column vector. 



P 7 o EXP ( i • PxRpd) (a)P2 



(97) 



(98) 



(99) 
(100) 



The above orbits are permutation equivalent to the 7x7 family found by 
Petrescu They all are cognate, and the P^ orbit is even self-cognate since 



P 



f) (a) = (pW(a))- 



(102) 



Also 

P 7 (1) (-a) = P T ■ P 7 (1) (a) • P for P = [e u e 4 , e 5) e 2 , e 3 , e 7 , e 6 ] (103) 

so P 7 (1) (-a) ~P 7 (1) (a) 

Due to some freedom in the construction of family components the method 
of Petrescu allows one to build other families of Hadamard matrices similar to 
P 7 ^. Not knowing if they are inequivalent we are not going to consider them 
here. 

5.7.3 The 'cyclic 7 roots' matrices 

There exist only four inequivalent 7x7 complex Hadamard matrices, inequiva- 
lent to P 7 and the 1-parameter family found by Petrescu (see previous subsec- 
tion), associated with nonclassical cyclic 7-roots. This result was obtained in 
[2*1] and is based on the catalogue of all cyclic 7-roots presented in |32 |. 
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The two matrices C jA , C\g correspond to the so-called 'index 2' Solu- 
tions to the cyclic 7-roots problem. They have the circulant structure [U] i ■ = 
mod 7)+i) where 

x = [1, 1, 1, d, 1, d, d] /or (104) 

x = [1, 1, 1, d, 1, 3, d] /or (105) 

and 

d= -3 + u/7 suchthat rf 2 + 3 d+1 = = ^ d .g =1 ( 106 ) 
4 2 

The corresponding dephased matrices are denoted as C 7 ° A and C 7 °g : 
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d 2 


d 


d 


There holds Cf^ 




r<(o) 

°7A 


( <= 


, rr(o) 

^ °7B — 


) and the matrices C 7A \ 









Òg} and C^j, C^j are equivalent to (c$) , (<5$) and fc§)) , ((^2) 



respectively (see l|17jl). 



The structure of CÍq, CÍf, related to 'index 3' solutions to the cyclic 7-roots 
problem is again [U] i . = x^_j moc j 7)+i, where x is a bit more complicatcd. 
We put: 

x = [1, A, B, C, C, B, Af for C$ (109) 
x= [1, A, B, C, C, B, A] T for C$ (110) 
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where A = a, B = ab, C = abc are products of algebraic numbers et, b, c 
of modulus equal to 1 (see the final remark 3.11 of 21 ). Their numerical 
approximations are given by 



a w exp(i ■ 4.312839) 
b w exp(i ■ 1.356228) 
c w exp(i • 1.900668) 



and then 



Again, C. 



.4 
B 
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a w (-0.389004) - 
ab w (0.817282) 4 
aòc w (0.280434) 



-i - (-0.921236) 
i • (-0.576238) 



i • (0.959873) 
denote the dephased versions of C. (0) 
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7( ., C*^' wn i cn read: 
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and 
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The matrices C 7C > and are symmetric and are related by complex con- 
jugation, CÍfy = C^q. All four cyclic 7-roots Hadamard matrices are isolated 

na. 



5.8 N = 8 

5.8.1 Orbits stemming from Fg 

The only maximal affine Hadamard family stemming from Fs is the 5-paramctcr 
or bit: 

F^\a, b, c, d, e) = F 8 o EXP(i ■ R„ m (a, b, c, d, e)) (121) 

where 
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R „(5)(a,ò, c, d,e) 
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(123) 



The family F$ is self-cognate. 

The above orbit is permutation equivalent to the orbit constructed with the 
Dinà's method 
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(124) 
(125) 



is the only maximal afiinc Hadamard family stcmming from Í4 and D(a3, a^, a§) 
is the 4x4 diagonal matrix diag(l, e la3 , e lai , e 1 " 5 ). 

Eq. (|124|) leads to F 2 ® F 4 for cti = . . . = 015 = 0, which is not equivalent 
to F 8 , see EHJ. However, F 8 (5) (0, 0, (1/8)2tt, (2/8)2tt, (3/8)2tt) is permutation 
equivalent to Fg, since 

F s = F 8 (5) (0, 0,(l/8)27r, (2/8)2tt, (3/8)2tt) • [e 1; e 3 , e 5 , e 7 , e 2 , e 4 , e 6 , e 8 ] T , (126) 

where ej denotes the z-th column vector of the Standard basis of C 8 . Thus 
Eq. (|124fl generates the only maximal afhne Hadamard family stemming from 
permuted F&. 

The matrix F^ (ir/2, n/2, 0, 0, 0) yields the only real 8x8 Hadamard matrix, 
up to permutations and multiplying rows and columns by — 1. It is dephased, 
so it is permutation equivalent to F 2 ® Fi ® F2 . 

Therefore all appropriately permuted tensor products of Fourier matrices, 
F<2 ® F2 <S> F2 , F2® F4 and F§ , although inequivalent [SÜ] , are connected by the 
orbit J32U. 
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5.9 iV = 9 

5.9.1 Orbits stemming from F 9 

The only maximal affine Hadamard family stemming from Fg is the 4-parameter 
or bit: 



where 
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with tu = exp(i • 2tt/9), so w 9 = 1, and 
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(129) 



The orbit Fg^ is self-cognate. Observe that its dimension is equal to the defect, 
d(Fg) = 4, which follows from Eq. ifTül) . 



It is permutation equivalent to the 4-dimensional orbit passing through a 
permuted Fg, constructed using the Dinà's method: 
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where D(a, (3) is the 3x3 diagonal matrix diag(l, e la , e 1 "). 

The matrix F 9 (4) (0, 0, 0, 0) = F 3 ® F 3 is not equivalent to F 9 [Sü|, but 



(130) 



F 9 = F 9 (4) ((l/9)27r, (2/9)2rr, (2/9)2tt, (4/9)2tt) • [ei, e 4 , e 7 , e 2 , e 5 , e 8 , e 3) e 6 , e 9 ] 



(131) 
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where ei are the Standard basis column vectors. 

Thus both inequivalent permuted matrices, F3 (8 F3 and Fg , are connected 
by the orbit fT5ü)l . 



5.10 N=10 

5.10.1 Orbits stemming from Fio 

The only maximal affinc Hadamard famílies stemming from Fio are: 
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(133) 
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(135) 



The affine Hadamard families F^q and (-^ío^ are cognatc. 

At least one of them must be permutation equivalent to an orbit constructed 
using the Dità's method, either: 
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where D(ai, . . . , CK4) is the 5x5 diagonal matrix diag(l, e 1Ql , . . . , e 1 " 4 ), or 

F$ B (0 1 ,...,0i) (137) 

such that its i, j-th 2x2 block is equal to [F^] i . ■D(a)·F<2, where i,j € {1 ... 5}, 
-D(a) is the diagonal matrix diag(l, e la ) and a = 0, /3i, . . . , /?4 for j = 1, 2, . . . , 5 
respectively. 

This is because F^ A (0) = F 2 ® F 5 and F{q' b (0) = F 5 ® F 2 are, according 
to |59|. permutation equivalent to -Fio, so both Ditjà's orbits are maximal affinc 
Hadamard families stemming from permuted -Fio's. 



5.11 jV = ll 

5.11.1 Orbits stemming from Fu 

The Fourier matrix Fu is an isolated 11 x 11 complex Hadamard matrix: 
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(138) 

where w = exp(i • 27r/ll), so w 11 = 1. 



5.11.2 'Cyclic 11— roots' matrices 

There are precisely two ~ equivalence classes of complex Hadamard matrices, 
inequivalent to Fu, associated with the so-called nonclassical 'index 2' cyclic 
11-roots. This result is drawn in |23 p.319. 

The classes are represented by the matrices C^ A , Ò[ 1B below, their re- 
spective dephased forms are denoted by C[° 1A , Cub- Both matrices have the 
circulant structure \U] i - = mo d u)+1j where 

x = [1, 1, e, 1, 1, 1, e, e, e, 1, e] for ò[^ A (139) 
x = [1, 1, e, 1, 1, 1, e, e, e, 1, e] for Cg^ (140) 
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and 
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Applying transposition to C^ A , cf^ B - equivalently to c[^ A , C[° 1B - yields 
a matrix equivalent to the original one (see l|47| l). 
The matrices are given by: 
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5.11.3 Nicoara's 11 x 11 complex Hadamard matrix 

Another equivalence class of 11 x 11 Hadamard matrices is represented by the 
matrix communicated to the authors by Nicoara: 
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where 




The above matrix is isolated 7 since its defect íí(A^ 1 ( 1 ) ) = 0. 



7 R. Nicoara, private communication 
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5.12 TV =12 

5.12.1 Orbits stemming from F í2 

The only maximal affine Hadamard famílies stemminj 

F 12 o EXP (i • 
F 12 o EXP (i • 
F12 o EXP (i • 
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= F 12 o EXP (i 
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where w = exp(i • 2n/12), so w 12 = 1, w 6 = — 1, = i, and 
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Thus we have three pairs of cognate famílies, and the F^ A family is self- 
cognate. 



At least one of the above orbits is permutation equivalent to the orbit con- 
structed using the Dinà's method: 

F 1 ( 2 9) (a 1 ,...,a 9 )= (158) 





[F 3 } 12 ■ £)(a 4 ,a5,a 6 ) • F[ 1 \a 2 ) 


[F 3 ] 13 • D(a7,a 8 ,ag) ■ Ff'(a 3 ) 


[F 3 ]^-F^( ai ) 


[F 3 ] 22 ■ D(a 4 ,a 5 ,a e ) ■ F 4 (1, (a 2 ) 


[ F 3] 2 , 3 • D(a r ,a s ,a 9 ) ■ F 4 w (a 3 ) 


. m 3A -F^( ai ) 


[F 3 ] 3 2 • D(a 4 ,a 5 ,a 6 ) ■ F 4 (1, (a 2 ) 


[ F 3} 3 , 3 • D(a 7 ,a 8 ,ag) ■ F 4 (1) (a 3 ) _ 



where D(a,j3, 7) denotes the 4x4 diagonal matrix diag(l, e 1Q , e 1/3 , e 17 ) and 
F 4 (1) (a) is given by ijH^I. 

It passes through F^ 2 \o) — F 3 tg>Fi, which is, according to [Sü], permutation 
equivalent to i*i 2 , so -F 12 is a maximal afhne Hadamard family stemming from 
a permuted F\2- 

It also passes through F 12 (n/2, n/2, 7r/2, 0), which is dephased, so it is a 
permuted F 3 ® F<z 8) -Fa • Thus permuted and inequivalent í\ 2 and i 7 ^ <8> F 2 (g) _F 2 
(see [5ÏÏ]) are connected by the orbit of (|158ll . 

Note also that a similar construction using the Dinà's method with the role 
of F 3 and F± exchanged yields a 7-dimensional orbit, also passing through a 
permuted _Fi 2 : F4 (g) F 3 , which is a suborbit of one of the existing 9-dimensional 
maximal affine Hadamard families stemming from F4 <g> F 3 . 



5.12.2 Other 12 x 12 orbits 

Other dephased 12 x 12 orbits can be obtained, using the Dità's method, from F 2 
and dephased 6x6 complex Hadamard matrices from section 15.61 for example: 
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where D(/3%, . . . denotes the 6x6 diagonal matrix diag(l, e 



ifi 



5.13 iV=13 

5.13.1 Orbits stemming from i*x3 

The Fourier matrix F13 is an isolated 13 x 13 complex Hadamard matrix: 
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where w = exp(i • 27r/13), so w 
5.13.2 Petrescu 13 x 13 orbit 



(2) 

There exists a continuous 2-parameter orbit P 13 of 13 x 13 complex Hadamard 
matrices found by Petrescu |62| . 



P$ (e, /) = P 13 o EXP (i ■ R p g (e, /) 



(169) 
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where í = exp(i • 2tt/30), so í 30 = 1, í 15 = -1, and 
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where 



G (/) = arg(-22^) + i. ^77-^(27)) - f (172) 

Since the function G(f) is nonlinear, the above family is not an afhne 
Hadamard family, but it is not clear whether it could be contained in any afhne 
Hadamard family of a larger dimension. 

Due to some freedom in the construction of family components the method 
of Petrescu allows one to build other similar families of Hadamard matrices. 
Not knowing whether they are inequivalent we are not going to consider them 
here. 

5.13.3 'Cyclic 13 roots' matrices 

There are precisely two ~ equivalence classes of complex Hadamard matrices, 
inequivalent to Í13, associated with the so-called 'index 2' cyclic 13-roots. This 
result is drawn in [21] p.319. 
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The classes are represented by the matrices C[^ A , C[® B below, their re- 
spective dephased forms are denoted by C^ A , c[^ B . Both matrices have the 
circulant structure [U] i ^ — x^j mo( j i3)+i, where 

x = [1, c, c, c, c, c, c, c, c, c, c, c, c] for C 1 -^ (173) 
x = [1, d, d, d, d, d, d, d, d, d, d, d, d] for c[^ B (174) 

and 



-iwn +1 .M5oT2^ (l75) 



1-^ , . fVl30-2^5\ 



V 12 / V 12 

Conjugating cjg]., fc = A, B yields a matrix equivalent to the original 
one PT) . 

The matrices C-[g^, and C^g, C^g^ are symmetric. They read 
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5.14 TV =14 

5.14.1 Orbits stemming from F14 

The only maximal affinc Hadamard famílies stemming from F14 are: 



F$(a,b, c, d,e,f) 
(F[f{a,b,c,d,e,f)) T 



F14 o EXP (i • R p (6) (a, b, c, d, e, /)) (181) 
• (R F m (a, b, c, d, e, /)) 



= F14 o EXP 



where 
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where w — exp(i • 27r/14), so w 14 = 1, w 7 = — 1, and 



(182) 



R m (a,b,c,d,e,f) = 
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F$ and (f^J are a pair of cognatc familics. 

At least one of them can be obtained by permuting one of Dinà's construc- 
tions, either 



p( 6 ) ( n 
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,ae) is the 7x7 diagonal matrix diag(l,e 1Q1 , . . . , e 1 " 6 ), 



where D (ai, 
or 

Fub(Pi, ■ ■ ■ ,/3e) (185) 
such that its i, j-th 2x2 block is equal to [F 7 ] i . •£)(«)• F 2 , where i, j € {1 ... 7}, 
-D(a) is the diagonal matrix diag(l, e la ) and a = 0,Pi, . . . , (3q for j = 1,2, ... ,7 
respectively. 

This is because F^(0) — F 2 ® F 7 and F^(0) = F 7 ® F 2 are, according 
to |59|. permutation equivalent to i 7, 14 , so both Dinà's orbits are maximal affine 
Hadamard families stemming from permuted -F'14's. 



5.14.2 Other 14 x 14 orbits 

Other dephased 14 x 14 orbits can be obtained, using the Dinà's method, from F 2 
and dephased 7x7 complex Hadamard matrices from section I5~7l for examplc: 
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where fc e {A, B, C, £>} and 

(i, m) e {(A, A), (A, B), (A, C), (A, D), (B, B), (B, C), (B, D), (C, C), (C, D), (D, D)} 



designate 7x7 complex Hadamard matrices associated with cyclic 7-roots, and 
D(fti,... ,(3%) denotes the 7x7 diagonal matrix diag(l, e i/?1 , . . . , e i/3e ). 



5.15 TV = 15 

5.15.1 Orbits stemming from F15 

The only maximal affine Hadamard famílies stemming from F15 are : 

F$ (o, b, c, d, e, /, g, h) = F 15 o EXP (i • R p w (o, ò, c, d, e, /, 3, fe)) (191) 
(F 1 ( f(o,6,c ) d ) e,/,5,ft)) T = F150EXP (V (i? F (8,( a ,ò,c,d, e ,/,.g») T ^ 
where 
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where w = exp(i • 27r/15), so w 15 — 1, and 



R F w(a,b, c, d,e,f,g,h) = 



The famílies and í-F'ig·'J are cognate. 
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At least one of them can be obtained by permuting one of Dinà's construc- 
tions, either 
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is the 5x5 diagonal matrix diag(l, 
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F$ B (/?!,..., /? 8 ) (195) 

such that its ï, j-tl·i 3x3 block is equal to [ísjjj • D(a,j3) ■ Fs, where i,j £ 
{1 . . . 5}, D(a, (i) is the diagonal matrix diag(l, e 1Q , e 1/3 ) and 

(a,/3) = (0,0), (/3!,/3 2 ), C9 a ,/3 4 ), (/%,#»), (/3 7 ,/? 8 ) 

for j = 1, 2, . . . , 5 respectively. 



This is because F$ A (0) = F 3 ® F 5 and F-^(O) = F 5 (g) F 3 are, according 
to |59|. permutation equivalent to -F15, so both Ditjà's orbits are maximal affinc 
Hadamard families stemming from permuted íis's. 
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5.16 iV = 16 

5.16.1 Orbits stemming from i<\ 6 

The only maximal affinc Hadamard family stemming from Fi§ is the 17-parameter 
orbit: 

F^ 7 \a,b, c, d,e,f,g,h,i,j,k,l,m,n,o,p,r) = F w o EXP(i • R ( i 7) (a, . .. ,r)) 

(196) 

where 
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where w = exp(i-27r/16), so u> 16 = 1, w 8 = —1, u> 4 = i, and (where typographic 
purposes we denoted e — a + k by 




so this is not a binomial coefficient) 
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The 17-dimensional family F^ 6 of JV = 16 complex Hadamard matrices is 
self-cognate. Note that this dimensionality coincides with the defect, d(Fie) = 
17, which follows from Eq. l(TÜ|) . 



The above orbit is permutation equivalent to the orbit constructed with the 
Dinà's method: 



" [F 2 ] ltl -FP( ai ,. 


• ,"5) 


mi,2 


D (au, ■ 


■ ,an) ■ F^ 5 \a e , . . 


■ ,a w ) 


[F 2 ] 2A -F^( ai ,. 


■,<*5) 


m 2 , 2 


D (au, . 


•>"i7) • Fs 5) ( a 6, ■ ■ 


■ ,aw) 



(199) 

(5) 

where the only maximal affine Hadamard family stcmming from F%: i*g such 

that Fg 5 ^(0) = Fg, is given by l|121[) . and D(au, • ■ ■ , «17) is the 8x8 diagonal 
matrix diag(l, e iai1 e iai7 ). 

This is because 

Fia = # 1 ( 6 17) (0,(l/16)2 7 r,(2/16)2Tr,...,(7/16)27r)· (200) 

[ei, e 3 , e 5 , e 7 , e 9 , en, ei 3 , ei 5 , e 2 , e 4 , e 6 , e 8 , ei , ei 2 , eu, e 16 ] T 

where denotes the i-th Standard basis column vector, so indeed (|199fl generates 
the only maximal affine Hadamard family stemming from a permuted F\e. 

~ ( 171 

The above orbit also passes through F{ 6 ; (0) = F 2 (8 -Fsj as well as through 
permuted i*2 ® -P2 (B> F4 and i*2 <8 F 2 <8> -F2 (8 -F2, since ^ 8 (5) of (PIjl . or the 
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permutation equivalent F§ ' of (|121|l . both pass through permuted F 2 ® F4 and 
F 2 (g> F 2 (g> F 2 . Note that all the tensor products iq 6 , F 2 (g>F 8 , F 2 (g>F 2 (g> F A and 
F 2 ® F 2 ® F 2 ® F 2 are inequivalent [21] • 



6 Closing remarks 

Let us summarize our work by proposing a set of dephased representatives Q n 
of equivalence classes of Hadamard matrices of size TV = 2,3,4,5, and by enu- 
merating the sets from the sum of which one should be able to extract such a 
set of representatives for TV = 6, . . . , 16. The dots indicate that the existence 
of other equivalence classes cannot be excluded. For instance, one could look 
for new inequivalent families for composite TV using the construction by DrÇà 
|22| with permuted some of the component families of Hadamard matrices of 
smaller size. 

We tend to believe that the compiled list is minimal in the sense that each 
family is necessary, since it contains at least some matrices not equivalent to all 
others. However, the presented orbits of Hadamard matrices may be (partially) 
equivalent, and equivalences may hold within families as well as between them. 

In the list of complex Hadamard matrices presented below, let {Xffi} denote 
the set of elements of the family xj^ . 



n = 2 g 2 = {F^'}. 

N = 3 g 3 = {F 3 (0) }. 

TV = 4 G 4 = {F 4 (1) (a);ae [0,tt)}. 

TV = 5 g 5 = {F 5 (0) }. 

tv = 6 g 6 C 
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U(l,m)e{(A,A),(A,B),(A,C),(A,D),(B,B),(B,C),(B,D),(C,C),{C,D),(D,D)}{^^lilm} 



N = 15 



Gu c 



{F™} U 



{(^) T } 



u 



N = 16 



016 C 




Note that the presented list of cquivalcnce classes is completo only for N = 
2, 3, 4, 5, while for N > 6 the full set of solutions remains unknown. The list of 
operi qüestions could be rather long, but let us mention here some most relevant. 

i). Check if there exist other inequivalent complex Hadamard matrices of size 



ii) . Find the ranges of parameters of the existing N = 6 families such that all 

cases included are not equivalent. 

iii) . Check whether there exists a continuous family of complex Hadamard 

matrices for N — 11. 

iv) . Investigate if all inequivalent real Hadamard matrices of size N = 16,20 

belong to continuous families or if some of them are isolated. 

v) . Find for which TV there exist continuous families of complex Hadamard ma- 

trices which are not affine, and which are not contained in afhne Hadamard 
families of a larger dimension. 

vi) . Find the dimensionalities of continuous orbits of inequivalent Hadamard 

matrices stemming from Fn if N is not a power of prime. 

Problems analogous to i)— iv) are obviously open for higher dimensions. 
Thus a lot of work is still required to get a full understanding of the prop- 
erties of the set of complex Hadamard matrices, even for one-digit dimensions. 
In spite of this fact we tend to believe that the above collection of matrices will 
be useful to tackle different physical problems, in particular these motivated by 
the theory of quantum information |43) . Interestingly, the dimension N = 6, 
the smallest product of two different primes, is the first case for which not all 
complex Hadamard matrices are known, as well as the simplest case for which 
the MUB problem remains open |54j . 
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